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Abstract 

We study vrvr scattering in the 1 = 1, = 1^ channel at finite tem- 
perature in the framework of the extended Nambu-Jona-Lasinio model 
that explicitly includes vector and axial-vector degrees of freedom in 
addition to the usual scalar and pseudoscalar sector. The S'-matrix 
in the coupled channels qq and vrvr is constructed via /3-exchange in 
the s-channel. The self-energy of the p-meson contains both quark 
and pion loop contributions. The analytic structure of the ^-matrix 
for T > is investigated and the motion of the p-pole as a func- 
tion of coupling constant and temperature is followed in the complex 
•\/i-plane. For numerical calculations, parameters are chosen in order 
that m-j^, /tt and the experimental vrvr phase shifts 61 at zero temper- 
ature are reproduced, and then the behavior of the p-pole as well as 
the vrvr cross section is investigated as a function of the temperature. 
We find that the position of the p mass stays practically constant for 
< T < 130 MeV, and then moves down in energy by about 200 MeV 
for 130 MeV< T < 230 MeV. 
25.75.-q, 12.39.Ba, 21.65.-hf, 24.85.-Fp 

Keywords: vrvr scattering, p-meson, finite temperature, analytic struc- 
ture. 
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I. INTRODUCTION 



During the last decade many experimental and theoretical efforts have been ded- 
icated to the study of relativistic heavy-ion collisions, in which one hopes to observe 
the two aspects of the QCD phase transition: (i) the confinement / deconfinement 
transition in which hadronic matter decomposes into the quark-gluon plasma, and 
the converse process, and (ii) that of chiral symmetry restoration / breaking, both of 
which should occur at high values of the temperature and / or density. Depending 
on the energy density reached, the hot fireball formed after the collisions may (or 
may not) enter into the quark-gluon plasma phase and undergo a chiral transition 
simultaneously before it cools down and hadronizes, leaving it a challenge for people 
to identify its trajectory from the observed hadronic products. In particular, during 
the evolution of the fireball, a large number of pions is produced in the mid-rapidity 
region, and one can expect that interactions between them play an essential role 
in the dynamical evolution. It is thus of interest to study the vrvr scattering in a 
thermal environment in order to understand the evolution dynamics of the system. 
The TTTT scattering cross section, for example, is an important input to transport 
equations that are used to simulate the dynamical evolution of the fireball. Based 
on hadronic models, Refs. have investigated the vrvr cross section in nuclear 
matter, and have claimed to observe considerable medium effects. In itself, a study 
of the temperature dependence of the vrvr scattering cross section is necessary, and 
we will examine this here in the framework of a chiral model. On the other hand, 
TTTT annihilation is one important source of dilepton production. In recent studies 
by the CERES and IIELIOS-3 collaborations 0, an enhancement of the dilepton 
spectrum has been observed in the energy range of 200 - 800 MeV in nucleus-nucleus 
collisions. One of the hypotheses brought forward to explain these data is that the 
p-meson mass shifts to lower values as a function of increasing temperature and 
density, and its width broadens while doing so 0]. Since a major source of dileptons 
in this energy range comes from the p-meson, it is necessary to understand to what 
extent a mass shift can occur, and if one occurs, whether it is an increase or de- 
crease. Thus, via our study of the vrvr system, we also wish to examine the medium 
dependence of the p-meson. These two problems, that of vrvr scattering at energies 
500 MeV< <900 MeV and that of the medium dependence of the p-meson, are 
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closely intertwined, and we address both of these problems in this paper. 

The medium effects of the p-meson have been extensively discussed recently and 
for a review we refer the reader to Ref. [Q. The p-meson properties at finite tem- 
perature are not yet well understood, and in particular the issue as to whether the 
p mass increases or decreases with increasing temperature has not yet been settled 
In Ref. 0, for example, a vector dominance model that does not support 
a chiral or deconfinement phase transition is used to describe the temperature de- 
pendence, and these authors find a weak rising temperature dependence of the p 
mass and width. On the other hand, recent lattice calculations have found no sign 
of a temperature dependence in the meson properties for temperatures up to 0.92Tc 
||10|| . In this work, we study the temperature dependence within a chiral model that 
displays a phase transition in the chiral limit and a smooth crossover for small values 
of the current quark mass. We investigate the analytic structure of the p propagator 
and consider the p-meson as a pole in the complex plane. In this way we contribute 
additional insight. 

Our starting point is an effective chiral model: we use a version of the Nambu- 



Jona-Lasinio (NJL) model [ p!T] , p!2| that has been extended to explicitly include vector 
and axial-vector degrees of freedom ||T3|^r5[ in order to study the vrvr cross section 
at finite temperature and vanishing chemical potential under inclusion of the p- 
meson. The NJL-model, constructed essentially on the observation of the chiral 
symmetry of quantum chromodynamics (QCD) gives a transparent description of 
the mechanism of spontaneous symmetry breaking. It has been extensively used to 
study the static meson properties, and has enjoyed remarkable success, especially 



in the pseudoscalar sector |[Tl|jr^ . However, to describe the p-meson in the context 
of the NJL-model is a delicate task |[T6|jr7| , since, on heuristic grounds, nip ^ 2mq, 
where ntg ~ 300-400 MeV is the mass of the constituent quark, and the NJL-model 
does not confine quarks. In the language of scattering theory, the p mass is close to 
the qq threshold. While theshold effects often modify the shape of a resonance on 
the real energy axis close by, the position of a pole in the complex plane is much less 
distorted. For this reason we choose to define the properties of the p-meson as those 
of the pole in the complex plane. The self-energy of the p-meson in the NJL-model 
arises from the quark loop contribution at the lowest order in an expansion in the 
inverse number of colors, l/Nc-, which leads to the unphysical decay p —>■ qq that is 
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due to the lack of confinement. Q On the other hand, it is known from experiment 
that almost all p- mesons decay via the channel p tttt, and thus a satisfactory 
description of the p-meson cannot be obtained unless the physical p-meson decay 
channel has been taken into account. From the point of view of the NJL model, the 
necessary pion loop corrections can be classified as part of the next order in the 1/Nc 
expansion |[T^j2^. In fact, these terms should represent the leading contributions 



at this order |^. It is also necessary to include the pion loop contribution to the 
p-meson self-energy in order to fulfill the unitarity of the tttt scattering amplitude. 
With the inclusion of both pion loop and quark loop contributions to the p-meson 
self-energy, the tttt scattering via p-meson exchange becomes a coupled channel {i.e. 
TTTT and qq channels) process. Care has to be taken to obey current conservation. 
We study the analytic structure of the p-meson propagator or scattering amplitude 
in the complex energy plane, in order to identify the p mass as a pole. This part of 
the work is completely new. The parameters of the model can be chosen such that 
at T = 0, the values of m^, and the experimental data of the tttt phase shifts 
are well reproduced. Temperature effects are incorporated via the imaginary time 
formalism and the temperature dependence of the tttt cross section is obtained. Note 
that TTTT scattering has been studied by several authors in the context of the NJL 
model, however only to obtain the scattering lengths at zero temperature without 
the inclusion of the vector meson sector and also with it |2^. Scattering lengths 



at finite temperature were calculated in p4| 



One may well ask the question as to the validity of the NJL model in describing 
the p-meson in the manner described. No conceptual difficulty arises with the cutoff. 
It restricts the three momentum |p| < A and these is no condition on the mass. The 
problem lies rather in the fact that the unphysical p —>■ qq channel is open, and the 
poles of the propagator, if taken only in the lowest order in the l/N^ expansion, will 



only reflect this physics. In spite of this, several authors p!8| , p!5[] have found that the 
static properties in the vector meson sector can still be well understood, and that 
this model can still be used in a range which may a priori seem to lie outside of 
its realm of validity. Although we have gone beyond this restrictive approximation 



^In spite of this deficiency, this model still appears to give a good description of the 
position of the p-mass [18,1^]. 
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allowing in addition the physical decay channel p —>■ mv to open, we are still able to 
observe and isolate the effects of the opening of the unphysical qq threshold. This 
can be clearly seen in the calculated vrvr cross section, and is discussed at that point. 

This paper is organized as follows. In Sec. 0, we describe our formalism in the 
context of an extended NJL-model for mesons as quark-antiquark excitations, and 
calculate the vector and pseudoscalar mesonic modes that correspond to the p-meson 
and TT-meson. We then calculate the pirn vertex which is an important element in 
the pion loop correction and also for tttt scattering. The physical process p —>■ it tt is 
then incorporated into the p-meson self-energy as a higher order correction in the 
inverse number of colors l/N^ in the NJL model. The calculations of tttt scattering 
amplitude and cross section are outlined in Sec. |T|. Numerical results are presented 
in Sec. where we show the variations of the position of the p-pole with coupling 
strength and temperature, and the temperature dependence of the vrvr cross section. 
In Sec. |V|, we summarize our results and their implications. Some calculations of 
quark loop contribution to the self-energy of different mesonic modes are collected 
in Appendix |^. 



II. MESONS IN THE EXTENDED NAMBU-JONA-LASINIO MODEL 

Our starting point is the extended two-flavor Nambu-Jona-Lasinio Lagrangian 
density: 



i^i'-T-i^f + i^iVi^f , (2.1) 



with the spinors i/) in Dirac, color and flavor space, and mg a small common current 
mass for the up and down quarks, tUq = ttIq = itiq. The r°, a = 1, 2, 3, represent the 
isospin Pauli matrices, and Gi and G2 are coupling constants of dimension [Mass]~^. 
Eq. is an extended version of the original NJL model which also includes 

terms associated with the vector and axial- vector meson modes. In the chiral limit, 
niQ = 0, this Lagrangian displays SUl{2)®SUr{2) chiral symmetry even for different 
coupling constants Gi and G2 in the scalar-pseudoscalar and vector-axial-vector 
sectors. The NJL model is not a renormalizable field theory due to the fact that 
the interaction between quarks is assumed to be point-like. It is therefore necessary 
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to specify a regularization scheme in order to define the NJL model completefy. 
Throughout this work, we shall employ a three-momentum cutoff scheme in which 
we introduce a parameter A, such that p < A. This cutoff A, together with tuq, Gi 
and G2, are the four parameters in our present model. Note that this Lagrangian 
has already been studied in detail elsewhere [IT^IIRO. We follow the notation of 



Ref. [^] and refer the reader to this paper for extensive detail. Therefore in the 
following subsections, we only present those details which we directly require for the 
calculations of the vrvr scattering in the p-meson channel. 

A. Mesons as quark- ant iquark excitations 

Mesons as qq bound states or resonances can be obtained from the quark- 
antiquark scattering amplitude. For example, in the present model, there are col- 
lective modes carrying the quantum numbers J^'" = 0^^, 0^^, 1 , and 1"^"*" that 
are composed of interacting qq pairs, and which have the same quantum numbers 
as the 71, a, p, and ai mesons. For brevity, we shall refer to these modes as meson 
states of the corresponding name in the following discussions. 

For our purpose of studying tttt scattering, we require the effective propagators 
of the TT and p-mesons, as well as the effective vertex of the nqq interaction, all of 
which can be found by constructing the quark-antiquark scattering amplitude. We 
list these results here first. 

Consider first the vector meson exchange channel, since in the present model, this 
channel is simplest as it does not couple with other channels. In the random phase 
approximation, the quark-antiquark scattering amplitude shown in Fig. |^ can be 
expressed to leading order in l/N^ as an infinite sum of quark- loop chains, and can 
also be recast into the form of a Schwinger-Dyson equation. The p-meson exchange 
amplitude D^^(g^) for quark-antiquark scattering can be generated by the one- loop 
polarization function nJ^J^^^(g^) via the Schwinger-Dyson equation of Fig. 0, 



[-tD^^n = h2tG,S^,gn + [-2zG,6^,g^'\ [-tU^^ [-zD^] , (2-2) 

where a, b, c, d are isospin indices, and p, u Lorentz indices. The notation VV refers 
to the vector-vector channel. The calculation of the one-loop polarization function 
n^J^^^(g^) is straightforward |15], and is collected for completeness in the Appendix 
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together with some other one-loop polarization functions for different channels that 
we need in the calculations below. The Lorentz and flavor structure for the p- 
meson channel allows for the decomposition into transverse (T) and longitudinal 
(L) components, i.e. 



n 



S. 



abi 



(2.3) 



where L^^ = q^cf and T^^ = g^^ — q^q^ with q^ = q^/ \f(^. Consequently the p- 
meson exchange amplitude can also be decomposed into transverse and longitudinal 
components as 



DT{q^)T^' ^Di^{q^)L 



where 



-iD 



-2iG2 



T,L\ 



l + 2G2U^Xiq^ 



(2.4) 



(2.5) 



It will be seen later that only the transverse exchange amplitude of the p-meson 
enters into the vrvr scattering amplitude. 

In the pion channel, the situation is more complicated because of the mixing 
between the pseudoscalar (vr) and axial- vector (ai) channels. In analogy to the case 
in the p-meson channel, the Lorentz and flavor structure of the relevant one-loop 
polarizations |TH] (see the Appendix) leads to the following decompositions. 



D^^.,{q') 



[D^\q')T^, + Dt^iq')L^,]5at 

D^^{q')q,Sat, 

D''\q^)q,5a,. 



abi 



(2.6) 
(2.7) 
(2.8) 
(2.9) 



where the indices A and P refer to the axial and pseudoscalar channels, respectively. 
The Schwinger-Dyson equations become a set of coupled equations that can be 
written in matrix form as 



iD = iK + iK [-in] [-iD] , 



(2.10) 



where 
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iD= , K= . 2.11 

This matrix equation has the solution 



-iD = iK(l-nK)-^ 
_ 2iGi ( A iB 
~ I -iB C 



(2.12) 



where A = 1 + 2^20^^, B = 2iG2UP^, C = -^2(1 - 2G'in^^)/G'i, and V = 
Det(l -nK). 

In order to resolve the mixing between pseudoscalar and axial- vector modes and 
identify the corresponding independent eigenstates and their coupling to the quarks, 
the matrix of the coupled exchange amplitudes which appears in the quark- antiquark 
scattering amplitude, 

- ^r ^ (.j.r. ^ _ .^^^.^ _ .^^^^^^ j ^ ^^^^^ j , (2.13) 

has to be diagonalized. One obtains the form 

-tT = -t%-t%, (2.14) 

with —i% and — iT^ being the qq scattering amplitudes for the pure pseudoscalar 
(tt) and axial- vector (oi) exchange, respectively. For our purpose, we examine only 
and approximate it near the pole position as 

- ~ [{ignqq " IQpy- )i75T] [{igTvqq + WpvTy )n'5^]^^ ^, (2.15) 

from which we can read off the complete Trgg vertex to be 

{ignqq " Wpv-;^)ilbT- (2.16) 
In this expression, the relevant coupling constants g.„qq and Qpy are given as 

(2.17) 



2 

9nqq 



and 
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gl, = 2G2(1 - 2GiU 



pp\ 



f2.18) 



The additional pseudovector coupling constant Qp^ is induced by the pseudoscalar- 
axial- vector mixing p^ . 

Note that in this section, our formulae have been explicitly stated for temperature 
T = 0. The generalization to finite temperature leaves these equations formally 
unaltered - in that case, however, Lorentz invariance is broken and the functions are 
no longer a function of g^, rather of and (f separately. The derivatives required 
in Eqs.( p.l7D and ( p.l8|) are then to be taken with respect to q^ and evaluated 
at q^ = m^, and = 0. However, in the following, for brevity we will still use 4- 
momentum argument in the expressions even for finite temperature, and they should 
be understood in the above sense. 



B. yOTTTT vertex 



The pvrvr vertex is an essential element which is required for constructing the pion 
loop graph and the vrvr scattering amplitude. As will be seen in the next subsection, 
the p-meson may couple to pions via a triangular quark loop which is shown in 
Fig. ^ Using Feymann rules for finite temperature, one may write 

-tT^;np',p) = 2 X (-1) X |^_0-Tr[(z^?.,, + z^?,,^)z75z5(A;) 

xiWnqg - Wpv^^)il5'i'S (p' + k)'jf,iS (p + k)] , (2.19) 

where the trace Tr[- ■ ■] is to be carried out over the color and spinor indices, f3 
is the inverse temperature, and the Matsubara sum in n runs over the fermionic 
frequencies Un = (2n + l)7r//9, with n = 0,±1,±2,---, which via k = (ic(j„,k) 
enter into the quark propagator S{k). Note that the pqq vertex is represented in 
Eq. ( p. 191) only in its spinor structure 7^. The coupling strength gpgg will always be 
included via the p-meson exchange amplitude into the scattering amplitude, which is 
our convention. In deriving Eq. ( |2.19| ), we point out that both graphs of Fig. |^ give 
equal contributions, since they differ from each other for the Dirac trace of the quark 
loop only by a minus sign, which will be cancelled out by another minus sign arising 
from the flavor factors. This is the origin of the factor 2 in Eq. ( |2.19| ). Note that 
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the flavor factors liave been omitted in Eq. (p.l9 ), and will be incorporated later. 
We have also introduced the three-momentum cut-off A to regularize the divergence 
of the integral. 

Direct evaluation of Eq. ( p.l9| ) and an analytic continuation to real variables 
leads us to the form 



- iVp:^{p\p) = -tV^^{p',p){p' + p) 



(2.20) 



for the vertex, where the dimensionless function V^tttt can be expressed in terms of 
several elementary integrals Ji, l2{p' — p), I^ip'jp), whose definitions are given in 
Eqs. (|A10a|) - (|A10c|) of the Appendix. One finds 



+ 



— Ami. 



(g2 - 2ml)[h{p' -p)- hip)] + 2m%{p',p) 



12ml 



{q^ + 2ml)h{p' -p)-2h 



(2.21) 



where the momentum of the p-meson is denoted by g = p — p', and is the number 
of colors. 



C. 1/Nc corrections 



In this section, we extend the calculation of the p-meson (vector-vector) polar- 
ization to the next order in the 1/Nc expansion. There are two types of diagrams 
that enter at this level, and they are shown as (b) and (c) of Fig. 0, i. e. we construct 

n*;r = n(fj + nw (2.22) 

where = nj;^, ^°°p) + nj^^ and li^^^ ^°°p) and nj^^ refer to the contributions 
of Fig. ^ (b) and (c), respectively. The function n|fj refers to the vector- vector 
polarization of a single quark loop as already introduced in Eq. (|2.3|) , Il^^j = IlJ^J^, 
and which is given in the Appendix, see Eqs.( |All| ) and ( [A12| ). The additional term 
that contains the pion loop is necessary on physical grounds, since the decay 
p — > TTvr is the main decay channel. Within the context of the NJL model, counting 
arguments in the expansion in powers oil/Nc allows one to classify this term as being 
of the next order in 1/Nc relative to the quark bubble. One notes that in addition. 
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the contribution n|j^ ^^'^^ is also of the same order. It is required on physical grounds 
to ensure current conservation [^]. Note that if one were to contract the quark 
loops of Fig. 11(c) to a point, one would recover the tadpole diagrams of Ref. We 
will also see in Sec. |T| that this choice of diagrams enables us to guarantee unitarity 
of the TTTT scattering amplitude. With the introduction of this extra contribution, 
the quark-antiquark scattering amplitude of Eq.( p.5| ) becomes modified to read 

~2iG2 



1 + 2G2 



(2.23) 



In the pole approximation to the pion propagator, the p self-energy arising from 
the pion loop as shown in Fig. |5(b), is given by: 



p-TTTTj 



X 



(2.24) 



Here the pion momentum p is in the imaginary time formalism p = (i/i„,p), with 
the bosonic Matsubara frequencies /i„ = 2mi/l3, n = 0, ±1, ±2, ■ ■ -. After evalu- 
ation of the Matsubara sum over n, the external p momentum q = {iUn, q) is to 



be analytically continued to construct a causal propagator, 



m: 



ie. Due 



to the complicated structure of the pirir vertex given in Eq. ( 2.201 ), however, the 
complete calculation of the real part of 11*^'^ is quite difficult. We shall instead 
approximate in Eq. ( p.24| ) the complete [—iT^{p',p)] by [—iVpT,T,{p,p)]{p + p')^ with 
{qo/2, ugQ/A — ml), and go the zeroth component of the p momentum. The 



p 

piTTT vertex is then independent of the integration variable and can be factored out 
of the integration. We find the transverse and longitudinal parts of the pion loop 
contribution to have the form 



(-g^ + 4m^)/W(g) + 2/l 



r(^) 



n 



(tt loop) ^^2 



-2i [-iVp^^ {p, p) fp^^] ^/i'"^ , 



(2.25) 
(2.26) 



where the integrals I^^^ and l2^\p) have the same structure as Ji and l2{p) given 
in Eqs. ( |A10a|) and (|A10b|) in the Appendix, except that the quark mass m is to be 
replaced by the pion mass m^, and the fermionic frequencies Un = {2n + l)7r//3 are 
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to be replaced by the bosonic frequencies fin = ^mr/jS with n = 0, ±1, ±2, ■ ■ -. An 
analytic expression for the imaginary part of 11^^°°'''' (g^) can also be derived. One 
finds 

Imni: '°°^\q') = -[-^V,UP,P)fp..?T^iqo - ^<f" coth(/3go/4), (2.27) 

where q is the momentum of the p-meson, which can be off-shell. 

One notices immediately that 11^ '""^""(g^) 7^ 0, as was the case for the quark 
loop contribution (see Eq. A12 ). A direct calculation of Figs, ^(c) for Ii'^^^^'^'> would 



ensure that 11^^ = |2^. We do not do this explicitly here, but evaluate this quan- 
tity indirectly assuming the constraint 11^ '""^""(g^) 7^ 0. Since II^^^,^^'' = U^'^ ^^'^^Ofiu, 
it follows that 11*^'^ ^'^'^^ = —11^ loop) rjj_^^g polarization functions are known. 



D. p-meson mass and analytic structure of the scattering amplitude 

In order to identify the mass of the p-meson, one is required to find the poles 
of the vector-vector scattering amplitude D'^'^ or the corresponding propagator for 
this channel in the complex energy plane. This translates to finding the poles of 
the transverse amplitude Dxiq"^) that is given in Eq.( p.23D . We thus investigate the 



analytic structure of Dt from Eq. ( p.23D in this section. 



To avoid the complexity associated with the analytical continuation of the pvrTr 
vertex in the complex plane, we shall assume a constant value for the dimensionless 
pnir vertex function V^^^r in the following study of the analytic structure, which is a 
simplification that will not change the feature of our results. 

The poles of the scattering amplitude are given by the zeros of the denominator 
of Dxiq"^), which is, according to Eq. ( [^.23| ), 

F{q) = 1 + 2G2 [n5?^(g') + ni^\q^)] . (2.28) 

Before we discuss the analytic structure of the function F{q), it is useful to first 
examine the quark loop contribution Ili^\q'^). The analytic structure of Il^\q'^) is 
governed by the integral /2(g), 

i d'^k 1 

ilM - . X -Ee-/ (^ |,._,^;| p + ,).-,„a - (2.29) 
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In back-to-back kinematics, il2 depends only on the zeroth component of the mo- 
mentum q. After evaluating of the Matsubara sum and performing an analytical 
continuation, one finds the integral in Eq. ( p.29| ) to be 



1 Mb JE^ - ml tanh((3E/2) 



where A^; = y^A^ -|- m^. 

We consider ihiQo) as a complex function in the complex go-plane. This function 
has branch points at go = ±2mg, ±2Ae, but for our purpose only the branch point 
at go = will be relevant. Choosing in the complex go-plane a cut starting from 
the branch point go = 2mq and running to the right, one has two regions on the first 
Riemann sheet in which il2{%) is to be evaluated separately: 

(i) On the first sheet of the complex plane away from the cut, where the integrand 
in Eq. (|2.3CI| ) is non-singular, ihi^o) is readily calculated as 



Ab JE^ - ml tanh(/3E/2) 



which can be expressed analytically at zero temperature by 



IT=0 



, A^ + A V4m2-g2 ^ 
in arctan 



(2.32) 



(ii) On the real axis for the values of go with 2mg < Re(go) < 2A£;, the integrand 
of ihiqo) develops a singularity at = go/2 if e ^ 0+. Making use of the formula 

lim = V-+ m6{y) (2.33) 



e^o+ y — le y 



one obtains 



1 Mb JE^ -mltanh(f3E/2) i Jx^-Amf, 

where V stands for the Cauchy principal value, and x the real part of go. At zero 
temperature the integral in Eq. ( p^.34D can be analytically evaluated as. 



ilo(x) 



lT=0 



167r2 



Ap-^A Jx^-Am^ xA- KeJx^ - Ami 

2 In — + ^ In ^, 

mq X xK + Ke^x"^ - Ami 



x^ — Am'l 



IQ-R X 



(2.35) 
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In order to study poles of the scattering amplitude in the complex plane which are 
located close to the physical sheet, one has to continue the analytic functions onto 
the unphysical sheet. The analytic continuation of il2{(lo) through the cut into the 
second sheet is obtained by p6| 



0) 1 2nd 



0) list 



4m2 - g2 



Stt 



■ tanh 



(2.36) 



where "1st" and "2nd" denote the first and the second sheet. This analytic contin- 
uation is sufficient for our investigation since we are only interested in the poles on 
the unphysical sheet that is close to the positive real axis of the first sheet. 

The integral l2^^ occurring in the pion loop contribution 11^^ (g^) (see Eq. (|2.25|) ) 
has the same analytic structure as the quark loop integral I2 discussed above, except 
that the quark mass and the hyperbolic tangent in I2 are replaced in Jg'^'' by the 
pion mass and the hyperbolic cotangent, respectively. 

The analytic structure of the function F{q) in Eq. (|2.28| ), whose zeros correspond 
to poles of the scattering amplitude, is determined by the integrals I2 in Il^\q'^) 
and Iq^^ in Il!j^\q'^) which are both two-sheeted functions. Consequently F{q), as 
a combination of II^?'* and \ has two branches into multi-sheeted functions. In 
Fig. H, we illustrate the structure of F{q) in the complex qo-plane: F{q) has two 
relevant branch points at go = 2m^, and at go = 2mg, and two branch cuts on 
the real axis starting from 2m^, and 2mq respectively and running to the right; A 
continuous path is drawn to show the joining of Riemann sheets along the branch 



cuts. We use the notation 



with i,j = 1,2 to denote different sheets, e.g. 



stands for the case that I2 in 11^'' takes value on the first sheet, and I2 



in n^''' on the second sheet, and so on. In order to familiarize the reader with the 
analytic structure. Fig. | shows a path through the multi-sheeted plane: starting 



from point A in the plane 



1(9) iM 



crossing the cut into the plane 



and 



following the path until at B one crosses again a cut into 



2(9) 



then turning 



around to C and entering into 



2(9) 2^'^^ 



etc.. 



The numerical evaluation of F and the ensuing determination of the pole position 
of DT{q^) for the p-meson mass, is discussed in Sec. IV. 



III. vrvr SCATTERING 
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A. s-channel p-meson exchange 



Within the extended NJL model, there are several types of diagrams that con- 
tribute to TT7C scattering. The full complement of processes can be found in Ref. ||23|| . 
Among these are the box and a meson exchange diagrams, which dominate the 
s-wave scattering and quantitatively reproduce the scattering lengths, but do not 
suffice to account for the full scattering amplitude as a function of the center of mass 
energy in the vector-isovector channel . 

In this work, our interest lies in particular in studying the physical situation in 
the mid-rapidity region in high energy heavy-ion collisions, where thermal pions are 
copiously produced with kinetic energy, i?^ ~ 200 MeV, which is of the order of 
the typical temperature scale reached by the system shortly after the collisions. At 
that stage, we expect tttt scattering to be characterized by a center-of-mass energy 
^/s > 2Et^ ~ 400 MeV. In this energy region, according to the phase shifts data of 



TTTT scattering |£8|, the vrvr cross section has already an important contribution from 
the p-meson exchange channel. Thus, in the following, we shall restrict ourselves to 
p meson exchange diagrams in the s-channel which contribute to the vrvr scattering 
at the lowest order in l/N^ in the extended NJL model as is shown in Fig. ^ The 
intermediate p-meson state in this graph is to be understood as being given by 
the full propagator corresponding to Fig. |^ and this is built up of quark-antiquark 
excitations as well as vrvr modes. Of all diagrams we evaluate only the s-channel 
contribution in Fig. ^ since it should be the dominant contribution to the cross 
section. The scattering amplitude for this graph is given as 

where pi, p2 and ps, are the momenta of incoming and outgoing pions, respectively, 
and /p^TT is the flavor factor to be associated with the pTrvr vertex. For example, 
ffm-n = 2 for a p^-K^-K^ vertex in Fig. H. At finite temperature, Lorentz covariance is 
broken, i.e., in general the cross section of a two-body reaction does not only depend 
on the temperature, but also on the cm. velocity of the initial pair with respect to 
the heat bath. In the following calculations we shall assume the initial pion system 
to be at rest in the heat bath. 

Inserting the Lorentz structure of the pTrvr vertex given in Eq. (|2.20| ), the scat- 
tering amplitude in Eq. (|3.1j ) can then be recast as 
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- iMTTTT^nTT = -i[-iVpjr7T{~P2,Pl)fp7Tn] " u) , (3.2) 

where s, t, and u are the standard Mandelstam variables, s = {pi + ^2)^? t = 
{Pi ~ PzY) and u = {pi — PiY ■ Note that the scattering amphtude is proportional 
to {t — u) = (s — 4m^)cos^ = (s — 4m^)Pi(cos 0), which is consistent with the 
expectation that p-exchange is a p-wave scattering. Only the transverse component 
of the p-exchange amplitude contributes to the scattering amplitude due to the 
specific Lorentz structure of the pun vertex and the fact that pions in the initial- 
and final-state are all on the mass shell. This is in accordance with the fact that 
the poles of the transverse part of the amplitude Dt determine the position of the 
p mass. 

The differential cross section follows as 

da{s,t) 1 

dt 167rs(s — Ami' 

and then the total cross section by 



\M r, (3.3) 



to 

ais) = I dt 



^da{s, t) 
~dt 



1 + /b(v^/2) , (3.4) 



where to = 0, ti = 4m^ — s, /^(a;) = l/[exp(/3x) — 1], and the factor [1 + fB{\fs/'2)'f 
has been included in order to incorporate the Bose-enhancement effects in the final- 
state phase space. 



B. Unitarity of the p-meson exchange amplitude 

Unitarity at zero temperature provides a strong constraint on any scattering cal- 
culations, while the situation for finite temperature is less clear. It is thus important 
to confirm that the choice of diagrams evaluated for the vrvr scattering amplitude 
lead to a unitary 5— matrix at T = 0. The vrvr scattering amplitude in our approx- 
imation is to be calculated from Eq. ( p. 21) . We note that if Dx^q^) is restricted to 
contain the quark loop contribution H^?'* alone, no unitarity condition can be fulfilled 
as n^?^ becomes imaginary only above the qq threshold ^ > 2mq. This threshold 
is an unphysical artifact of the NJL model, at least in the confined phase. Including 
the term U^^ opens the vrvr threshold for ^/s > 2m^, and allows the optical theorem 
to hold for vrvr scattering. We recall once again that, from a physical point of view, 
this channel is essential, since nearly 100% of the p-mesons decay into pions. 
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It is apparent from this discussion that tttt rnr scattering within this model 
forms part of a coupled channels system that also involves the qq sector. The 
scattering amplitude for vrvr tttt proceeding via p-meson exchange is only one 
element in the two-channel S'-matrix, 



S 



Q _ - Q - 

q _ Q 

^TTTV — *qq ^nn — >7Tn 



(3.5) 



The other channels are all to be understood as proceeding via an s-channel p-meson 
exchange, and can be envisaged from Fig. 5(a) by simply constructing the relevant 
subprocess from this diagram. The matrix elements for the qq —>■ qq and qq ^ tttt 



are given as 



-iM 



qq-^qq 



-i[-iVpqqfpqqfDT{s){t - u) , 



and 



-iM 



qq—>-7V7r 



-i[-iVpqqfpqq][-iVp^^fpj,n]DT{s)(t - u) 



(3.6) 



(3.7) 



Thus, noting the explicit form of Eqs.( p.2|) , (|3.6| ) and ( |3.71 ) with Dt = 2G2/F, and 
F as given in Eq.( 2.28| ), one may factorize each term in the S matrix as 



I tAqqA^jF 

lA^^Aqq/F l + iAlJF 
In this expression, At^tt follows directly from Eq.( ^.2| ) to be 



(3. 



^2 ^ G2is - Ami) I s - Ami , ... . ,2 



12tx V s 
while similar arguments can be given to extract as 

Gois-Aml) I 



(3.9) 



A 



Am'i 



QQ 



127r V s 
Unitarity of the 5* matrix, SS~^ = 1, implies 

Al^ + A\^ = -2ImF = -AG^lm 



PQQi 



n. 



in) 



- AG2lm 



n 



(9) 



(3.10) 



(3.11) 



and similar result for Im 



n 



(<?) 



Making use of the explicit result for Im II^^^ from Eqs. (|2.27|) of the previous section 



one can actually show that A, 



2 

TTTT 



-4GoIm 



(tt) 
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and 



-4^21111 



n 



(9) 



T 



. Therefore, after including the pion loop contribution to 
the p-meson exchange amplitude, unitarity can be guaranteed. 

For later use, we recall that the elastic cross section in the l^^ partial wave 
scattering has a unitarity bound: 

4n(2l + 1) 



Pc 



(3.12) 



where Pcm is the cm. momentum. Since the p-meson exchange corresponds to a 
p-wave (/ = 1) scattering process, the cross section should always be bounded by 



-'cm' 



IV. NUMERICAL RESULTS 



In this section, we present our numerical results. We first discuss our choice of 
parameters, which are set at T = 0. We then discuss the numerical evaluation of 
the p-meson mass in the complex plane, and the tttt scattering cross section for a 
constant value of the pun vertex. A complete calculation, including the momentum 
dependence of the vertex function is presented in the ensuing subsection. 



A. Parameters 



Before discussing our numerical results, we have to specify the four parameters 
in the present model, namely the current quark mass mo, coupling strengths Gi, and 
G2, and the three-momentum cut-off A. Our choice of parameters uiq — 3.7 MeV, 
A = 650 MeV, GiA^ = 2.5, G2A^ = 2.3, leads to a value of the pion mass as 

= 140 MeV, pion decay constant /^r = 85 MeV, and quark condensate = 
—(267 MeV)^, which corresponds to a quark mass = 458.2 MeV at temperature 
T = 0. In particular, we have chosen the value of G2A^ in such a way that at 
T = 0, the p-meson appears as a tttt resonance at rup ~ 768 MeV. The dynamically 
generated quark mass nig resulting from this parameter set is deliberately chosen 
to be large, so that the qq threshold lies above the physical p-meson position and 
the p ^ qq channel is not open at T = 0. The value for the condensate (t/^-^) is, 
however, still acceptable. 
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B. /9-poles and cross section for constant \V, 



pTVTV I 



As in Sec. |II D| , a constant (energy-independent) value of the pmr vertex function 
\VpT,T,\ is used to avoid the intricacies induced by the comphcated structure of the 
pTTTc vertex. We have checked this assumption via a direct calculation whose result 
is shown in Fig. ^. For a rather wide energy range below the qq threshold, IV^ttttI 
depends only weakly on the energy. In this section, we will assume a constant value 
of the vertex, IV^ttttI = 2.4. 

Fig. 1^ shows the energy dependence of different real and imaginary parts of 
F{^/s), the denominator of the transverse p-meson exchange amplitude Dt{s) given 
by Eq. ( p.28| ), at the temperature T = 0. From this figure, the necessity of including 
the pion loop contribution 11^ into the p-meson self-energy becomes apparent. If 
there is no pion loop contribution to F{y/s), below the qq threshold ^/s = 2mg, the 
imaginary part of the quark loop 11^^ vanishes, since the decay channel p ^ qq is 
closed, and one finds a real root of -F(y^) at ~ 620 MeV where l-|-2G'2Re IIj?^ = 
0. This real p-pole of Dx{\/s) corresponds to a p-meson as a qq bound state, and 
would lead to an unphysical divergence in the vrvr cross section. 

If we take the pion loop contribution 11^ into account, the imaginary part of 
the polarization, Im II^'' -|- 11^^ ^ below the qq threshold, as long as it is above 
the TTTT threshold y/s = 2m^. At the position ^/s ^ 800 MeV, where the real part of 
F{^/s) vanishes, the imaginary part of F{^/s) does not vanish, and we have to move 
into the complex energy plane to find a complex root rrip of F{q), which corresponds 
to a p-meson as a vrvr resonance with Re(mp) ~ 800 MeV and finite width. 

Fig. II illustrates p-pole trajectories for different temperatures as a function of the 
coupling strength G2^'^- At T=0, for a very large (unphysical) value G'2A^ = 6.4, the 
p-meson appears as a bound state at 200 MeV and vanishing imaginary part since it 
is located below the vrvr threshold. With decreasing values of G2^^ the p becomes less 
bound, until above the vrvr threshold the trajectory moves away from the real axis 
and Im(mp) becomes negative and of the order of 100 MeV. A cross on the trajectory 
marks the physical value of 6*2 A^ = 2.3. Decreasing the value of G2A^ even more the 
trajectory approaches the real axis again, and enters into the sheet 2^'^^ l'^'^^ closely 
above the qq threshold. This exercise repeated for several temperatures is shown 
in Fig. 0. If one follows the position of the crosses, which indicate the physical 
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location of the p as a function of temperature, one reads off the main result of 
our investigation: very little changes for Re(mp) for temperatures between and 
130 MeV. This can be more clearly seen from Fig. where we show the temperature 
dependence of Re(mp), Re(m7r) and nig. The pion Mott temperature Tm, which is 
defined by mT,{TM) = '^rrigiTM), is found to be Tm — 280 MeV. This somewhat large 
Tm is a consequence of the large rrig ~ 460 MeV at T = 0. 

Fig. shows the cross sections for vrvr scattering, which are calculated for a 
constant IV^ttttI = 2.4. At T = 0, we find that the p-pole is located at Re(mp) ^ 
815 MeV and Im(mp) ~ —86 MeV. This p-pole corresponds to the peak of the cross 
section at T = in Fig. |10| with a width F ~ — 2Im(mp) ~ 170 MeV. As expected, 
the cross section at T = touches the unitarity bound. In comparison with the 
cross section at T = 0, the peak of the cross section at T = 150 MeV becomes lower, 
but its position does not shift much. While the cross sections at T = and 150 MeV 
have resonance-like shapes and can be directly related to the corresponding p-poles 



given in Fig. |10|, the relationship between the p-pole position and the shape of the 
cross section at T = 200 MeV is less clear. In the latter case, the p-pole on the 
li<i)^2^'^^ sheet is found to be slightly above the qq threshold that is indicated in 
Fig. ^. One may expect in this case, the shape of the cross section is a combination 
of the effects of the qq threshold and the p-pole we found. On the other hand, the 
cross section at T = 220 MeV displays mainly the threshold behavior, known in the 
literature as the Wigner cusp. In this case the p-pole shown in Fig. M locates on the 



2(9) 



sheet, and is so far away from the physical sheet that it cannot make its 



existence dominant in the shape of the cross section. 



C. Cross sections for energy dependent \Vr 



piTTT I 



In our calculations, we have chosen the parameter G2 so that at T = the 
p-meson appears as a vrvr resonance with mass rrip ~ 768 MeV, as one can see 
from a fit of the cross section cr(7r+7r^ tt+tt^) given in Fig. ^ for an energy 
dependent IV^ttttI- The data points of a{n^n^ tt+tt") translated from the phase 



shifts of vr+TT scattering in the vector-isovector channel |2^, instead of the tt+tt 



total elastic cross section {e.g., given in Ref. [^), has been used here in order 
to single out the dominant p-meson exchange contribution to cr(7r+7r^ tt+tt^) 
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and remove background contributions resulting from other channels that are not 



present in our calculations. One sees in Fig. |Tl] that the cross section is properly 
bounded by the unitarity at T = 0, and that both the position and the shape can be 
reasonably reproduced by the present parameter set. Fig. |Tl] also shows the results 
for a{TT^7c^ TT+vr^) for a constant value of the vertex ( |V^,r7r| = 2.4) and for an 
energy dependent vertex |Vp^7r|- The more realistic choice of the energy dependent 
vertex, with the exception of the area surrounding the qq threshold, gives a better fit 
to the experimental cross section, which has been calculated from the phase shifts 6} 
for vrvr scattering. The kink on the right hand side of the curve for energy dependent 
IV^ttttI is due to the opening of the unphysical p qq channel. 

Fig. |12| illustrates the tt+tt^ elastic cross sections at different temperatures for 
energy dependent IV^ttttI, which is entirely determined in Eq. (|2.21| ) by the present 
model. These results show analogous features to those for a constant |Vp^7r| (Fig. |10]) 
in that the position of the maximum moves to smaller values of y/s with increasing 
temperature. One notes that the height of the maximum decreases with respect to 



that shown in Fig. |T0| for a constant vertex. This is because in this case the cross 
section is always dominated by qq threshold effects, displaying Wigner cusps. 



V. RESULTS AND CONCLUSIONS 



To summarize this work, we have studied the p-meson propagator and vrvr scat- 
tering via s-channel p-meson exchange in the framework of the extended NJL-model. 
In particular, we have argued that in addition to the quark loop contribution, the 
p-meson self-energy should include a pion loop contribution from physical consider- 
ations which is found to be essential in order to preserve unitarity of the ^-matrix 
and to reproduce the experimental vrvr phase shifts at zero temperature. In order to 
maintain current conservation, other diagrams of the same order in 1/Nc expansion 
must be included. These diagrams reduce to the well-known tadpole diagrams of 
the hadronic model ||^ in the limit in which the quark loops are contracted to a 
point. 

The new part of our work is a careful study of the analytic structure of the 
two-channel scattering amplitude. We have studied the p-pole trajectories in the 
twofold cut complex energy plane with respect to the variations of the coupling 
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strength and temperature. In order to reproduce the data of tttt phase shifts at zero 
temperature, it is required that the unphysical qq threshold should be put above the 
physical p mass, and the p-meson can be consequently identified as a vrvr resonance 
below the qq threshold. At temperatures T < 215 MeV, the p-meson appears as a 
complex pole close to the physical sheet. Specifically, as the temperature increases, 
for temperatures T < 130 MeV the p-meson mass increases only slightly and its 
width is almost unchanged, while for temperatures 215 > T > 130 MeV, both 
the mass and width of the p-meson decrease. For these temperatures the vrvr cross 
section has a resonance- like shape, and its temperature dependence is rather weak 
and is mainly determined by the temperature-dependent p-meson mass and width. 
At higher temperatures (T > 215 MeV), however, the poles locate on the unphysical 
sheet and qq threshold effects determine the shape of the cross section. The rather 
weak dependence of the position and width of the p-pole as a function temperature 
observed in this work is in agreement with the results of lattice calculations |]ID| . 
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APPENDIX A: ONE-LOOP POLARIZATION FUNCTIONS 

In the imaginary time formalism for finite temperature, one loop polarization 
functions for different meson modes are given by 



-^K^M = ^E^^^""/ (^(-l)Tr[(7.rJ^5(A:)(7.r,)z5(A: + g)], (Al) 
= ^Ee^^"V ^(-l)Tr[(^75r.)^5(A:)(^75r.)z^(A: + g)], (A2) 
-^K» = ^Ee^^""/ 7|^(-l)Tr[(.75r.)^^(A:)(7,75r.).^(fc + g)], (A3) 
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j fA fj'~> U 

-<a.(?) = ^T.^'^'''J J^si-m[il,l^ra)zS{k){z^,n)zS{k + q)], (A4) 
-^KU<i) - ^12^'^"' j (^(-l)Tr[(7.75ra)i5(A;)(7.75r.)i-S(^ + g)]. (A5) 

Here Tr[- ■ ■] is a trace over the color, flavor and spinor space, f3 is the inverse 
temperature, and the Matsubara sums on n run over the fermionic frequencies 
Un = {2n + l)7r//3, with n = 0, ±1, ±2, ■ ■ ■, which through k = {iUn, k) enter the 
quark propagator S'(A;) = {f^ + mq)/{P -rrig) = [(icc;„)7o -k • 7 + mq)/[(ia;„)^ 
with El = k^ + m^. After performance of the Matsubara summation the Umit rj ^ 
is to be taken, and an analytic continuation is to be applied to the boson momen- 
tum q — {i/im, q), A*m = 2m7r//9, m = 0, ±1, ±2, • • •, which implies ifim —>■ qo- The 
three-momentum cut-off A is used to regularize the divergences of the integrals. 

These polarization functions can be decomposed according to their flavor and 
spinor structure as, 

KU^")= l^ViQKu + nr^(g) V] Sa,, (A6) 
KUl")- ['^T^il")T,. + ni%')L,^]Sa,, (A7) 
nir,(?^) = n^^(g^)g,5„,, q, = (A8) 

n„r(g^) = n^^(g^)<5„„ (A9) 

where the standard longitudinal and transverse tensors are L^,^ = and T^,^ = 
9nu - q^qu- 

In terms of the following elementary integrals 



i d^k 1 

(27r)3 - m^] [{k+pf-mf ^^^^"^^ 
i d^k 1 

^3(p',p) = ^E^"""; (27r)3 (F - ml) [{k + p')' - "^^] [(^ + P? - ' ^^^^""^ 

one-loop polarization functions can be evaluated as 

nr(?') = ^^^^/A^c[(?' + '2ml)h{q) - 2hi (All) 

nr(g') = 0, (A12) 

n^^(g2) = AtNfN.h - 2tq^NfNj2{q), (A13) 
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n^^(g2) = Am,NfNj2{q)/^Jq^ (A14) 

n^^(g2) = -4m,NfNMq)/^, (A15) 

nf^(g') = -SimlNfNMq), (A16) 

n^^(g2) = nr(g^) + n^^(a. (ai7) 
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FIGURES 




FIG. 1. Quark-antiquark scattering in the random phase approximation, and its re- 
cast in the context of the Schwinger-Dyson equation. Thick wavy hnes indicate the full 
propagator of /9-meson, and thin wavy lines the bare propagator. 




FIG. 2. pTTTT vertices to the lowest order in l/Nc- the two diagrams differ by the 
direction in the quark loop. They contribute equally as explained in the text. 
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TT 




FIG. 3. Quark loop and pion loop contributions to the p propagator: (a) represents 
the leading order contribution in l/N^ expansion, while (b) and (c) are the next order 
corrections. Thick wavy lines indicate the full propagator of p-meson, and thin wavy lines 
the bare propagator. 
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FIG. 4. Joining of Riemann sheets of -F(go) along two branch cuts on the real axis 
starting from thresholds 2m,r, 2mq and running to the right. Different sheets are denoted 



by 



with i,j = 1,2 representing the first, second sheet in qq and tttt channels, 



respectively, (see also text.) 



28 



TT 



Pi 



P3 



7T 



IT 



P2 




Pa 



IT 



FIG. 5. TTTT scattering via s-channel p-meson exchange to the lowest order in 1/Nc 
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^ (MeV) 

FIG. 6. Energy dependence of the pvr-Tr vertex given by Eq. ( 2.21| ) at temperature 
T = 0. The discontinuity at -y/i ~ 920 MeV relates to the qq threshold 2mq. 
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FIG. 7. Energy dependence of real and imaginary parts of the denomenator of the 
p-meson propagator F in Eq. (^.28|) at T = 0. The breaks at ^/s ~ 920 MeV are due to 
the qq threshold. 
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FIG. 8. Variations of the p-pole position in the complex energy plane with the coupling 
strength G2A^ at different temperatures. The trajectory at T = begins with = 6.4 

for Re(mp)=200 MeV, and it moves to the right with decreasing ■ Similar for other 
temperatures. For each temperature the qq threshold 2mq is indicated by a dotted-triangle, 
and the p-pole position for the physical parameter G2A^ = 2.3 is marked by a cross. 
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FIG. 9. Temperature dependence of Re(mp), Re(m7r), and rriq. The long-dashed hne 
indicates the values of Re(mp) corresponding to p-poles in the sheet 2^'^), l^'^) 
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FIG. 10. Calculated tttt cross sections in the vector-isovector channel at different tem- 
peratures for a constant vertex |Vp7r7r| = 2.4. The unitary bound 127r/p^^ at T = is also 
shown. 
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FIG. 11. Energy dependence of the vrvr elastic cross section in the p exchange channel at 
T = for two choices of the /jvrvr vertex function (full energy dependence and approximated 



by a constant), compared to the data (represented by crosses) taken from Ref. ||2^. Also 
shown is the unitary bound 127r/p^^ at T = 0. 
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FIG. 12. Elastic cr(7r+7r" vr+vr") at temperatures T = 0, 150, 180, 200, and 220 MeV 
for energy dependent |V^7r7r|. The crosses represent the experimental cross section at T = 0, 
and the long dashed line is the unitarity bound 127r/p^^ at T = 0. 
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